Fully Self-Consistent Hartree-Fock Calculation of Jellium Slabs: Exact Treatment of Exchange Operator

We present Hartree-Fock surface energies, work functions and dipole barriers for a Jellium slab model at different electron densities and slab widths. The fully self-consistent calculations take into account the nonlocal exchange coupling of the momentum parallel to the surface with the perpendicular component of the orbitals. Typical oszillations due to quantum size effects have been observed. Our results provide lower and upper bounds for the Hartree-Fock surface energy and work function of the semiinfinite jellium model which can serve as benchmarks for previously reported variational calculations performed with different types of non-self-consistent orbitals. Furthermore we observed the expected characteristic band structures and electron densities, i.e., vanishing density of states at the Fermi energy and pronounced Friedel oscillations at jellium surfaces. In order to resolve the perpendicular components of the orbitals we have used a wavelet based interpolating scaling function basis. PACS numbers: 71.10.Ca, 73.20.At, 31.15.Ne, 02.70.Dh ∗Present address: Institut für Mathematik, Technische Universität Berlin, Straße des 17. Juni 136, D-10623 Berlin, Germany